Josephson junctions in narrow thin-film strips 
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We study the field dependence of the maximum current I m (H) in narrow edge- type thin- film 
Josephson junctions. The junction extends across thin-film strip of width IV C A = 2A 2 /d, the 
London depth A 3> d, d is the film thickness. We calculate I m (H) within nonlocal Josephson 
electrodynamics, which takes into account the stray fields affecting tunneling currents. For W <C 
ccf>o/8it 2 Ag c , g c is the critical sheet current density, the phase difference along the junction depends 
only on the junction geometry and the applied field, but is independent of the g c , i.e., it is universal. 
Zeros of I m (H) are equidistant only in large fields (unlike the case of junctions with bulk banks); 
they are spaced by ~ (j>o/W 2 that is much smaller than cf>o/W\ of bulk junctions. The maxima of 
I m (H) decrease as 1/yH, slower than 1/H for the bulk. 
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The physics of the edge-type thin-film Josephson junc- 
tions (e.g., two films in the (x,y) plane touching only 
along the edges at x = with no overlap) differs from 
that of the junctions with bulk banks mainly because 
of the stray fields, that affect the currents in the junc- 
tion and in the thin- film banks. The phase difference tp 
across the junction is also affected by the stray fields. As 
a result, ip is described by an integral equation, i.e., the 
problem becomes nonlocal [J 0, H,H, H, @] • 

Development of nonlocal electrodynamics of such junc- 
tions is still in progress and is a subject of grow- 
ing interest 0]. Long-range stray- fields are relevant 
for physics of sequences of interchanging 0- and ir- 
junctions [1, i, El EE 0, E, El- These anoma- 
lous chains of tunnel Josephson junctions are stud- 
ied also for the thin-film superconductor-ferromagnet- 
superconductor heterostructures 111, Il2fl. asymmet- 
ric grain boundaries in YI^CusO? y [lOL [l3|, and 
YBa2Cu307_ x /Nb zigzag junctions [14j . 

The phase distribution ip(y) along thin- film edge- type 
junctions has a length scale I = c0Q/87r 2 Ag c , the thin- 
film analog of the Josephson length [fj; g c is the critical 
sheet current density, A = 2X 2 /d, X 3> d is the London 
penetration depth, and d is the film thickness. We show 
in this work that when the width W of the junction con- 
taining strip (and the junction length that are the same) 
is less than £, the distribution of the phase difference <p(y) 
becomes £ independent, i.e., the same for junctions with 
different critical currents. In other words, for W -C £, 
f(y) is a universal function, that depends only on the 
applied field and the junction geometry. 

In this situation, we evaluate the field dependence of 
the maximum supercurrent I m (H) through the junction 
that turns out quite different from the standard Fraun- 
hofer pattern of bulk junctions. Zeros of I m (H) become 
equidistant only in large fields unlike in bulk junctions, 
and are separated by AH ~ <po/W 2 , which is much 



smaller than 4>o /WX of bulk junctions of the same length. 
The maxima of I m (H) decrease as 1/yH, that is signif- 
icantly slower than 1/H for the bulk. We show that 
Im (H) for a SQUID made of narrow thin-film strips with 
edge-type Josephson junctions differs remarkably from 
the canonic pattern of the bulk junctions. 

Let the x axis be along the strip and z be perpendicular 
to the film; the junction is located at x = 0, < y < W. 
The sheet current density g = (g x >9y) can always be 
written as g = curl Sz = (d y S , —d x S), where the S(x, y) 
is the stream function Since the current component 
normal to the edges is zero, S is constant along the edges 
(y = 0, W) and the total current through the strip is 



1= / dyg x = / dyd v S(0,y) = S(W)-S(0). (1) 
Jo Jo 

The London equation integrated over the film thickness 
reads: 

2ttA 



r-\\- 



■ C uvl z g=^S(x) i p'(y), (2) 

where h z consists of the applied field H and the part 
related to g by the Biot-Savart integral. The right-hand 
side here is a manifestation of a general rule: the field 
of a Josephson junction is equivalent to the field of a set 
of vortices distributed along the junction with the line 
density <p'(y)/2ir 

In narrow strips, the self-field due to the current g is of 
the order g/c, whereas the second term on the left-hand 
side of Eq. ((2]) is of the order gA/cW ^> g/c. Hence, the 
self-field can be disregarded, unlike the applied field H. 
Substituting curl z g = —V 2 S in Eq. |(5J), one obtains: 

2vrA 



^S^-^5(x) l p'(y) + H. 



c 2vr 
This linear equation has solutions S = S± 



(3) 
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So such that 



(4) 
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C(j)Q 



V 2 5 2 = — 



H 



(5) 



The boundary condition fl} is satisfied if Si(W) 
5i(0) = and S 2 (W) - 5 2 (0) = I. Hence we have: 



5i(r) 



dpS(u)^-G(r,p), 

Z7T 



(6) 
(7) 



Here, r = (x,y) and p = (u,v); G(r,p) is the Green's 
function for Eq. Q with zero boundary conditions that 
satisfies (2TrA/c<fio) V 2 G = — 6(r — p), an equation well 
studied in electrostatics 11611: 



G(r,p) = tanh 



sin V sin Y 



cosh(AT — U) — cos Y cos V 



(8) 



Q = 4tt 2 AG/c4>o, the capitals stand for corresponding co- 
ordinates in units of W/tt. The Green's function G(r, p) 
gives in fact the current distribution of a single vortex at 
r = p. 

Clearly, Si describes the current perturbation due to 
the junction. The first term in S 2 represents the screen- 
ing currents due to the applied field, whereas the second 
is due to the field of a uniform transport current. 

Given the stream function, we obtain the sheet current 
density through the junction: 



g c smip(y) = g x (0,y) = d y S(0,y) ■■ 

' W dv m dyG{0 ,y,0, v) ' U 



(9) 
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We rewrite this integral equation for the phase f{y) as: 



W 
where 



dVtp'jV) sinV 
cos Y — cos V 



h Y - 



8tt 2 A 3c ' 



8ir 2 A 

C(j} 



(10) 



(11) 



are the characteristic length, the reduced field, and the 
reduced current. 

To establish the boundary conditions for ip(y) we em- 
ploy the London equation for ^(±0, y) on the two junc- 
tion banks 



,(±o,y) 



~4tt 2 A 



dx(±0,y) 2tt 

— At. 



dy 



(12) 



where y) is the phase and A is the vector potential. 
We subtract these equations and utilize the continuity of 
A to obtain (p'(y) oc g y (0, y). The current g y must vanish 
at the junction edges, i.e., 



(f'(0) = <p'(W) = 0. 



(13) 
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FIG. 1: (a) The function tp'oifi) calculated according to 
Eg. (|16|l . (b) The solid line is <£>o(a0 obtained by numerical 
integration of tp'o(fi) shown in the panel (a). The dashed line 
is the approximation tpo(n) = 0.43 sin/i. 



We note that the length £ along with g c , the only 
material parameter of the junction, enters only the left- 
hand side of Eq. (flQ|) . In narrow junctions with W <C I, 
this term can be disregarded. While neglecting the term 
oc W/£ we have to disregard also the transport current 
i; otherwise, integrating both sides of Eq. (fTOf over the 
strip does not produce identity. 

The truncated Eq. (fTT)|) reveals a remarkable feature of 
junctions in narrow strips: the phase is just proportional 
to the applied field and can be written as ip(y) = h(po(y) 
where ifo{y) is an universal function governed by an in- 
tegral equation 



ip' (V) sinV 



dV 

o cos Y — cos V 



7T 

Y -2 



0. 



(14) 



which does not contain g c , the physical parameter of the 
junction quality. To study this function, we introduce 
s = cosl/, t = cosY and write Eq. (fl4|) in the form: 



1 



1 J(s)ds 
2vr 7_i t-s 
■ dip 



B n (t) . 



J = 271-^1 - s 2 -p , B n = - sin -1 1 . 
as 



(15) 



The reason for this manipulation is this: Eq. (|15p is the 
Biot-Savart expression for the normal component of the 
"field" B n at the surface of a thin strip — 1 < s < 1 car- 
rying the "sheet current" J(s). This integral equation 
can be inverted [171 ] . One, however, should have in mind 
that the current J(s) is not determined uniquely by one 
field component; currents of the form C/y/1 — s 2 with an 
arbitrary constant C correspond to full Meissner screen- 
ing and to zero normal component of the "field" . The 
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latter flexibility allows us to obtain the solution <p' (V) of 
Eq. (|14p that satisfies the boundary conditions (fT3]) : 



f 



7T Z COS /i 



zr/2 



T) cos 2 ?7 d?7 
-tt/2 sin// - sin 77 



(16) 



where the origin is shifted to the strip middle for conve- 
nience, fi = Y — 7r/2. 

The integral in Eq. (|16p is understood as Cauchy prin- 
cipal value and can be done numerically. The univer- 
sal function <p' {n) so calculated is shown in Fig. [T](a). 
The result of the numerical integration of this function 
obtained requiring (po((J,) to be an odd function of /1 is 
shown in Fig. [T](b). In particular, this calculation gives 
^o(tt/2) - W>(-f/2) « 0.86. 

Thus, for any applied field in narrow thin-film junc- 
tions the phase takes the form = hipQ^fi) + 0, 
where 6 is a constant. The total current through the 
junction is 



9cW 



r/2 



tt/2 



(17) 



Maximizing this with respect to provides 8 = tt/2 and 
the maximum current Z m : 



1 



g c W 7T 



tt/2 

(i/i cos[/i^?o(m)] 

-tt/2 



(18) 



Hence, I m (H) can be evaluated numerically; a good ap- 
proximation for I m (H) can be obtained as follows: 

The odd function fo(p) can be written as the Fourier 
series ^ a n sin(2n + l)/i to satisfy the boundary condi- 
tions (|13p . We take the lowest approximant cpo — ao sin /1 
with a = 0.43 to fit the difference <p (W) - <p Q (0) = 0.86 
that is found integrating numerically the exact derivative 
in Eq. (|16|) . The comparison of the phase found numeri- 
cally with ao sin /1 is shown in Fig. [T](b). 

In this approximation we have: 



9cW 



r/2 



d/i cos(/i do sin /1) 



tt/2 



= \J (a h)\. (19) 



Figure [2] shows that this approximation is quite accu- 
rate as compared to I m {H) calculated numerically with 
the help of Eq. (fT%|) . Zeros of the Bessel function Jo (a;) 
are equidistant for large arguments, but they are spaced 
roughly by it everywhere. Hence zeros of I m (h) are sep- 
arated by aoAh ~ tt, or in common units by: 



AH ~ 1. 



W 2 



(20) 



It is worth recalling that in bulk junctions of the length 
W the zeros are separated by AH w 2cj)o/W\ that ex- 
ceeds by much the thin-film spacing. 




FIG. 2: The maximum supercurrent i m = I m /gcW versus the 
normalized applied field h n — 4aoW 2 H/tyc^o. The dashed line 
is the approximation (1191) . 
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FIG. 3: Sketch of a rectangular SQUID made of two narrow 
thin-film strips with identical edge-type junctions 1 and 2. 



In the high-field region one can use the large argument 
asymptotics of Jq to obtain: 



0.61 g c 



cos 1.72- 



HW 2 



(21) 



Thus, the maxima of I m {H) decrease as 1/vH, i.e., 
slower than in the bulk case where I m oc 1/H. 

It is worth noting that in high fields the maxima I m (H) 
do not depend on the junction length W. Qualita- 
tively, this comes about because the tunneling current 
9x = 9c sinQvpo + 9) oscillates fast for h ^> 1 so that 
most of the junction length does not contribute to the 
total current, unlike currents in narrow belts of the width 
5 ~ O.3^/0o/ff near the strip edges. 

Let us consider now current flowing through rectangu- 
lar SQUID (superconducting quantum interference de- 
vice) made of narrow thin-film strips with two identi- 
cal Josephson junctions sketched in Fig. 03 In zero field 
the current distribution is symmetric with respect to the 
SQUID center and the line integral of g along any sym- 
metric contour is zero. When the field is applied, this 
symmetry is violated by the screening currents. How- 
ever, at the contour in the strips middle (shown in the 
figure) the screening currents vanish so that the contour 
integral of g remains zero. This contour crosses the junc- 
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FIG. 4: The maximum supercurrent i m — I m /2g c W versus 
flux <f)l4>a for a rectangular SQUID (Fig.© with Ao/W 2 = 5. 



tions at their middle, where the local coordinates fi = 0. 
Clearly, the flux (j> enclosed by this contour does not 
change if the contour is shifted as a whole by fx. In- 
tegrating g = —(c(j)o/4:ir 2 A)(Vx + 2TrA/(f) ) over such a 
contour we obtain: 

<P2<J*)-'Pi(j*)=tor4-' ( 22 ) 
The total current through the system is given by: 



e?^i(sin tpi + sin <p 2 ) 
2 

sin(hipo + 9) + sin ( h<po 



9cW 

rir/2 
dfl 

'-3T/2 

„ r /2 , . f, „ *<i>\ 

2 1 d/i sin nipo + 9 H cos — 

J-tt/2 \ 90 J \<P0 



2irc 



(23) 



As above, 9 is a constant with respect to which the cur- 
rent should be maximized. The maximum current then 
corresponds to 9 = ir/2 — ttq 



I m = 2g c W 



W 



J(\ 4an — ; — COS 7T — 



.4 



o (po 



(24) 



where Aq is the area of the "central" contour. Note that 
our argument is valid if the SQUID hole area is large rela- 
tive to the area of superconducting branches. In this case 
the difference between the flux enclosed by the "central" 
contour and the SQUID hole area can be disregarded. 

Thus, the standard SQUID pattern given by 
| cos(7r</>/</>o)| is modulated in our case by a slow vary- 
ing Bessel function. An example of I m (<j>/(f>o) is shown 
in Fig.Hfor a rectangular SQUID with A /W 2 = 5. We 
stress again that the pattern shown is obtained for large 
area SQUIDs made of narrow thin-film branches; for re- 
duced areas the interference patterns become more com- 
plex, a subject for further study. 



Summarising, we have evaluated the field depen- 
dence of the maximum supercurrent in narrow edge-type 
Josephson junctions in thin-film strips; the strip width 
W is supposed to be less than both the Pearl length A 
and the thin-film Josephson length £ of Eq. (fTTj) . Calcu- 
lations are done in the framework of nonlocal Josephson 
electrodynamics. We demonstrate that the stray fields 
cause a pattern I rn (H) with much reduced distance be- 
tween zeros, AH ~ cf>o/W 2 , and with a slow decreasing 
maxima in high fields, I m (H) oc l/y/H. The flux depen- 
dence of the maximum supercurrent through a SQUID 
made of narrow thin-film strips with edge-type junctions 
differs by much from the standard periodicity. 
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